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We investigate the statistical properties of wave functions in chaotic nanostructures with spin-orbit 
coupling (SOC), focussing in particular on spatial correlations of eigenfunctions. Numerical results 
from a microscopic model are compared with results from random matrix theory in the crossover 
from the gaussian orthogonal to the gaussian symplectic ensembles (with increasing SOC); one- and 
two-point distribution functions were computed to understand the properties of eigenfunctions in 
this crossover. It is found that correlations of wave function amplitudes are suppressed with SOC; 
nevertheless, eigenfunction correlations play a more important role in the two-point distribution 
function(s), compared to the case with vanishing SOC. Experimental consequences of our results 
are discussed. 
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I. INTRODUCTION 

Spin-orbit coupling (SOC) has the potential to make 
novel electronics applications possible^ as it allows one to 
control the electron's spin degree of freedom through its 
motion. Most systems of interest for such applications 
are nano- or mesoscopic in size, including semiconduc- 
tor quantum dots^ 2 - metallic nanoparticles^ and quantum 
corrals defined on surfaces A The energy spectrum, and 
more generally, properties of these systems are (typically) 
described by random matrix theory (RMT)^ 

In RMT, the system's properties are a consequence 
of its symmetries — in the classic Wigner-Dyson en- 
sembles, the key symmetries are time- reversal (T) and 
spin-rotation (er) invariance.«^ With both T- and a- 
invariance, the system is described by the gaussian or- 
thogonal ensemble (GOE); SOC breaks the cr-invariancc 
(while preserving T-invariance) , driving the system to the 
gaussian symplectic ensemble (GSE). [Systems with bro- 
ken T-invariance are described by the gaussian unitary 
ensemble (GUE).] More specifically, systems with SOC 
are described by random NxN matrices (with TV— ^oo) 
having quaternion components 



H = S®I 2 




(1) 



j'=i 



where S is an NxN symmetric matrix, and the {Aj} 
are NxN antisymmetric matrices; {c J } are the Pauli 
matrices, and I2 is the 2x2 identity matrix. A in Eq. [T] 
is related to the SOC of the microscopic Hamiltonian — 
A=0 in the GOE, while \=y/4N in the GSE. 

As most nanoscale systems of interest are de- 
scribed by RMT; 2 - - — it is important to understand the 
regimes/behaviors which arise with SOC and the prop- 
erties in these regimes. In this work, we consider the 
spatial properties of wave functions in (two-dimensional) 
chaotic nanoscale systems with SOC. The spatial prop- 
erties of wave functions often determine the system's re- 
sponse to experimental probes, and are important for 
devices/applications.— ~— While other works have dis- 



cussed properties/consequences of eigenvector statistics 
with SOC ) 15 ' 16 here we consider the spatial properties 
of eigenvectors and, in particular, how these properties 
evolve with the SOC. 

In what follows, we consider the properties of the 
Hamiltonian 



H = — p 2 + az-(pxa) + V(r) , (2) 
2m 



where V(r) is a confining and/or disorder potential. 
Results obtained via RMT arc compared with those 
obtained by direct simulation of Eq. [2] for a stadium 
billiard.— To characterize the system and understand 
its properties, one- and two-point distribution functions 
were computed in the crossover from the GOE to the 
GSE (with increasing SOC). In particular, it is found 
that excellent agreement between RMT and microscopic 
simulations are obtained in a "mean-field" description 
of the (GOE-GSE) crossover (see below). A key observa- 
tion from our results is that correlations of wave function 
amplitudes are suppressed with SOC. Interestingly, how- 
ever, these correlations play a more important role in the 
two-point distribution function(s), compared to the GOE 
(with vanishing SOC). 

The rest of the paper is organized as follows. The 
description of wave function statistics in RMT and, in 
particular, the description of the GOE-GSE crossover is 
discussed in Sec. [TTJ Details of our microscopic calcula- 
tions — namely the stadium billiard considered as well 
as the numerical approach employed — are presented in 
Sec. IIIII Our results are presented in Sec. IIVI — one- and 
two-point distribution functions obtained via RMT are 
compared with numerical results from the stadium bil- 
liard. Finally, Sec. [V] contains a summary of our results 
as well as remarks on experimental consequences. 
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II. WAVE FUNCTION STATISTICS IN RMT 

In RMT, wave function correlations are governed by 
the functional probability distributio n 18 ' 19 



= Afexp 



E 



drdr>:(r)tf s , s ,(r,r')Vv(r' 



(3) 

K s .s' (r, r') is the functional inverse of the two-point cor- 
relation function (tp* (r)tp s , (r')) , where the angular brack- 
ets (• • • ) denote an average with respect to the pa- 
rameter (3 depends on the system's symmetries — /3=1 
{)3=2) in the GOE (GUE), while /3=4 in the GSE. [N is 
a normalization constant.] V(ip) is the probability that 
a particular energy cigcnfunction with spin-er is equal to 
the specified function Vv(r). 

A key property of Eqs.[T]and[2]is their invariance under 
time-reversal; as a result, the energy levels are two-fold 
degenerate — the eigenstates {ip(r), Tip(r)} are degener- 
ate, where T is the time-reversal operator. Explicitly, 




(c) X 


*=0.09 


i J 






i k 


ft 











(d) 



X =0.25 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Y Y 



FIG. 1: (color online) Distribution of the {7^}, P({7i}) (from 
Eq. O, where A*=A v / 4iV. (a) A* = 0.05 (b) A* = 0.07 (c) 
A* = 0.09 (d) A* = 0.25. Notice all V change rapidly for 
A* < 0.1 and become sharply peaked at 7 ~ 1/2 for large A*. 
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As a consequence of this two-fold degeneracy, the wave 
function amplitude probed numerically and experimen- 
tally is |V'o-(r)| 2 = |(/>(r)| 2 -l-|x(r)| 2 - As noted above, we 
are interested in the regimes/behaviors which arise with 
SOC — we will not only be interested in the GSE, but 
also in the crossover from the GOE to the GSE. As such, 
we decompose the complex wave functions (f>(r) and x( r ) 
in Eq. 2] into their real and imaginary parts. Then, the 
wave function amplitude is parameterized as 

\Mr)\ 2 = ll 2 (r)+ 7 I ^(r)+7f xl(*)+ll xl(r), (5) 

where the parameters {7^}, which satisfy the constraint 
7i + 72 + T3 + it = 1) characterize the crossover — 71 =1 
with 7^=0 for i 7^ 1 in the GOE, while 7i=l/2 (i = 
1 • • • 4) in the GSE; in the crossover, the {7^} fluctuate 
and, hence, physical quantities must be averaged over 
their distribution. 

We obtained Vdji}), the distribution of the {7^}, nu- 
merically from Eq. [T] by considering the various orthog- 
onal invariants^ — the results are shown i n Fi g. [TJ We 
see that the V({"fi}) change rapidly for A/V AN < 0.1 — 
in particular, the P({7i}) are broad for small A, but be- 
come sharply peaked gaussian-like for larger values of A, 
moving towards 7^ = 1/2 (V i) with increasing A. Figure 
[2] shows the variance of the {7;}, (li)-(li) 2 , as a func- 
tion of A; the inset shows how the average values of the 
{ji}, (7s)) evolve with A. We see that the variance is ex- 
tremely small for larger values of A; even for small values 
of A (where the 'P({7i}) are broad and asymmetric), the 
variance does not exceed 0.03. As noted above, physical 
quantities must be averaged over the V({ji}); however, 
as will be seen below, rather good results are obtained in 
a "mean-field" description (due to the small variances), 
similar to what has been observed in the GOE-GUE 
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FIG. 2: (color online) Variance of {7^} vs. A*=A\/ 4N. Inset: 
average values of {7^} vs. A*. 



crossove r 10 ' 21 — rather good results are obtained by ap- 
proximating the {ji} by their average values (rather than 
averaging over the 7 3 ({7i})). 

From Eq. [31 all spatial correlations can be obtained 
once the two-point correlation function (tp*(r)ip s ,(r')) is 
known. To determine this, we expand the wave function 
as 



V>0) = ^2i'+A r ) c +.P +^-, P ( r )c-, 



p • 



(0) 



where the two-component spinors if) +tP (r) and ^+, p (r) 
are eigenstates of Eq. [2] with V(r) = 0. Explicitly, 
the eigenvalues {E + ,E-} and corresponding eigenstates 
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{^+,p(i-),V- pM} ar e (& = 1) 

£± = j^-iald; V± P (r) = --L= ( _ L . 1 /I . ^ e ipr 
2m r±,pv ^ ± tz /\z\ J 

where z=p x +ip y . The spectrum above describes two 
spin-split chiral surfaces with energy E, shown schemat- 
ically in Fig. [3K where k± = \/2mE + m 2 a 2 =F ma. 

To compute (tp*(r)i() s ,(r')), the Fourier coefficients (in 
Eq. |6|) are taken to be gaussian random variables having 
zero mean and variance given by 8 - (a, 6=+,—) 

« P c f,,k) = S a ,bSpMj^d(e{p) - e) , (c a , p Cf,, k ) = . 

(7) 

Writing the wave function as per Eq. [4] and using the 
parameterization in Eq. [31 we obtain^ (i,j=l,2) 

(4, i (r)4> j (r')) = {xi(r)x j (r'))=5 i , j f , (8a) 

(<Mr) Xj (r')> = -<Xi(r)^ (r')> = <^ 3 , (8b) 

where 

/ = ^[J (k+R) + Jo(k-R)}, (9a) 

g = ^[Ji{k + R)-J x {k-R)}. (9b) 

In Eqs. [9al and I9bl Jq(x) (J±(x)) is the Bcsscl function 
of order-0 (order-1)^ R=\r — r'|, and k± are the wave 
vectors associated with the chiral branches at energy E. 
The physics of Eq. [7] (and Eqs. [5a] and [SbJ) is that the 
system ergodically samples the energy surfaces 2 ^ (shown 
schematically in Fig. [3Ji). 

III. NUMERICS 

As described above, we are interested in comparing 
results obtained via RMT with those obtained by direct 
simulation of Eq. O To this end, we have computed the 
local density of states (LDOS) for a stadium billiard^ 
where the billiard's wall was constructed with a unitary 
delta-function potential^ 

V(r) = V 5(r-R(s)) , (10) 

with R(s) parameterizing the wall (and Vq — > oo). The 
retarded Green's function (GF) for the system, 

G(r,r';u J ) = (r\(u J -H + l 0+y 1 \v'), (11) 

is computed from the Dyson equation 

G(r,r';w) = G (r,r'; W ) 

+ V Q J^ds G (r,R(*);w) G(R(s),r';w) . 

In this equation, Go(r, r';w) is the free-particle GF, i.e. 
the GF in the absence of the corral's wall, but in the 
presence of SOC^. 

G (r,r» = G° (R;uj) I + G 1 (R;u>) ( .% ~ ie ^ ) 
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FIG. 3: (color online) (a) Spin-split energy surfaces with wave 
vectors k+ and (b) Stadium billiard considered in this 
work, (c) Spatial scan of the LDOS of a typical chaotic eigen- 
function. 

where 

G° (R;uj) = - t ^{k-H { 1 \Rk„) + k + H { 1 \Rk + )} , 

Gl(R;u) = -^{k^H{ 1 \Rk-)-k + H[ 1 \Rk + )} , 

(12) 

and exp(i#) = [(x — x') + i(y — y')]/R, with Hq 1 \x) and 

h[ 1 \x) being Hankel functions^ with R=\r — r'| and 
k± defined as before. The LDOS is then obtained from 
the GF via A(r,ui) = -(l/V)ImTr [G(r, r; w)]. 

The stadium billiard we consider is shown schemati- 
cally in Fig. [5b. With energy in units of Eo=l/(2mR l ) 
and SOC in units of a.Q=l/{mRo), where Rq is the ra- 
dius of the stadium's circular cap, we have considered 
cigenstatcs with energy i?~405£'o, and have investigated 
SOCs in the range 0<a< lOao- [Choosing i? =70A, and 
m=0.26m e (with m e being the electron's rest mass), one 
obtains ao=3.7 x 10 1 eVm, a value consistent with e.g. 
electrons on an Au(lll) surface ! 27 ' 28 ] A spatial scan of 
the LDOS for a typical eigenstate considered is shown in 
Fig. [3b; from the LDOS, one- and two-point distribution 
functions were computed, going from the GOE to the 
GSE (with increasing a). 

IV. RESULTS 

We now analyze the properties of the system, com- 
paring results from RMT with those obtained by direct 
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simulation of Eq. [2] for a stadium billiard. We begin by 
determining the regimes which arise as function of the 
SOC strength. To this end, we consider the one-point 
function Viv) = (S (y — A\ip a (r)\ 2 )} , which is obtained 
from Eq. [3] by integrating out the degrees of freedom ex- 
cept at r. Using Eq. we obtain 



V(y) 



v 



x exp 



x I 



471727374 
v 
4 
/ 1 



/' 

Jo 

(1-z) 



dz 



(13) 
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where Iq(x) is the modified Bessel function 
order-zero.— This expression reduces to "Pgoe(^) = 
exp(-v/2)/V2w in the GOE (71=1 and {7i}=0 for i^l) 
and Pgse(^) = 42/exp(-2i/) in the GSE (7, = l/2 Vi). 

Figure H] shows numerical results for V(v) for different 
values of the SOC; the results are compared with Eq.[T3l 
in a "mean- field" description, i.e. with the {7^} evaluated 
at their average values — for a=0.2ao (a=0.5ao), we 
find A = 0.04 VM (A = 0.08VlN)M The physics of 
Eq. [2] is determined by its two length scales — the spin- 
flip length / s [=l/(ma) and the linear dimension of the 
system L (~ Rq). Figure @] shows how the system evolves 
toward the GSE as the SOC is increased. In particular, 
we find the system to be in the GSE for a > 1.5ao i.e. 
l s { < 2i? /3; once the system is in this GSE regime, the 
statistics do not change further as the SOC is increased. 

We now turn to spatial correlations of cigenfunctions. 
We first consider the amplitude correlator C aa i (r, r') = 
(j4| (r) | 2 (r') 1 2 ) . Using the parameterization in 
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FIG. 4: (color online) V{v) = (6 (u - A\ip a {r)\ 2 )) . From top 
to bottom: a=0 (GOE), a=0.2, a=0.5, q=1.5 (GSE). Each 
curve has been vertically offset by one unit for clarity. 
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FIG. 5: (color online) Amplitude correlator C CTCT <(r, r') = 
(^|7/) CT (r)| 2 yl|V CT '(r')| 2 ). From top to bottom: a=0 (GOE), 
a=0.2a u , a=0.5ao, a=1.5ao (GSE). Each curve has been ver- 
tically offset by 1/4 unit for clarity. 



Eq. O we obtain 
<w(r,r') 



2[(7i 4 + 7 2 4 + 7 3 4 + 7 4 4 )/ 2 

+z{ihl + ihl)g 2 ] ■ 



(14) 



Notice that this reduces to C^° E (r, r') = 1 + 2/ 2 in the 
GOE, and to C^ s , E (r,r') = 1 + V 2 /2 in the GSE where 
V 2 =f 2 +g 2 . Numerical results for C CTCr /(r,r') are shown 
in Fig. [5] and are compared with Eq. 1141 with the {7^} 
evaluated at their average values (as before). We see that 
the maximum is larger in the GOE; more generally, the 
correlations decay more rapidly with SOC — amplitude 
correlations are suppressed as c-invariance is broken. 

Having determined the parameter regimes and, in 
particular, how large the SOC must be to be in the 
GSE, we now consider in greater detail the prop- 
erties of the system in the GSE. To this end, we 
consider the joint distribution function V(y\,v-i) = 
(<y(i<i-A|Va(r)| 2 )<y(i^-A|Va'(r / )| 2 )> 1 which is ob- 
tained from Eq. [3] by integrating out the degrees of free- 
dom except those at r and r'. For the GSE we obtain 



PgseK V2) = v 8 ( f^ 2) ex P (- 2 *sr) J i ( 4 <*» . ( 15 ) 

where Ii(x) is the modified Bcsscl function of order- lj^ 
for comparison, we also consider V(y\, V2) in the GOEi£ 

exp (-As/2) cosh (A» 

In the above equations, Xg={v\ + v 2 )/(l — X 2 ) and 
X P =X^yij{\ - X 2 ) where X=V (X=f) for the GSE 
(GOE). 

We now consider the properties and consequences of 
"P{ v ii v 2)- We begin by considering the conditional prob- 
ability 



(16) 



5 




which describes the wave function distribution at r 2 , pro- 
vided Al^ri)! 2 = u\. It follows from Eg. fl6l that cor- 
relations between fluctuations at different points depend 
on their amplitudes^ — regions of high amplitude (i.e. 
large v\) are correlated over larger distances, while re- 
gions of small amplitude are correlated over shorter dis- 
tances. P° SE (^ 2 ) for the GSE is shown in Fig. [5Ji for 
several values of V = ^Jf 2 + g 2 ; T^{v 2 ) for the GOE 
is shown in Fig. |BJd for comparison, for several values of 

/• 

From Eq. 1161 one can obtain the average (1*2)1/1 an d 
the mean squared fluctuation ((5^ 2 ) 2 ) ! / 1 = {v 2 ) Vi — (^2)^, 
where (■ • • ) Vl denotes an average with respect to V Ul (^2): 

(1/2)1/1 = 1 + ^2(^-1), (17) 
{{5v 2 ) 2 ) Ul = C[l + 2X 2 {v l -l) + X\l-2v l )} , 

where C=2 for the GOEji 2 . while C=l/2 for the GSE. [As 
before, X=V (/) for the GSE (GOE).] From this, we sec 
that fluctuations are suppressed in the GSE compared to 
the GOE. More generally, fluctuations are largest in the 
GOE (compared with the GUE^ and the GSE, Eq. [T7J 
and, hence, correlations are the weakest. 

We now consider the distribution of the product 
A|W(r)Vv(r')l, V(T) = (S (T - ,4|W(r)W(r')|))- ^( r ) 
determines a number of experimentally relevant quanti- 
ties, such as the form factor in resonant scattering in com- 
plex nuclciifii amplitudes in tunneling measurements, 
and the conductance amplitude distribution through 
small quantum dots.— From Eq. I15[ we obtain for the 
GSE 



|V|(1- V 2 



1-V 2 



1- V 2 



(18) 

where Kq(x) is a modified Bessel functions of order- 
zero^ in the GOE, we obtain 



^GOE(r) 



-Kr 



I-/ 2 



cosh 



/r 



Figure [7] shows results for ViT) for several values of 
V (/) for the GSE (GOE). We see that the maximum of 



FIG. 7: (color online) Product distribution V(Y) = 
(S(F - A|Vv(r)Vv(r')D) in the GSE for several values of V. 
For comparison, "pGOE(r) is also shown for / = 0.5. Inset: 
Comparison of numerical and RMT results for R = 0.055i?o- 



■PcsE(r) decreases with increasing V with the tail becom- 
ing slightly longer. For comparison, Vgoe(T) is shown 
for different values of /. We see that correlations play 
a more significant role in the GSE — indeed, except for 
a very small region near r=0, 'Pgoe(F) is essentially in- 
distinguishable from the result with /—►(). This is a con- 
sequence of the fact that fluctuations are largest in the 
GOE and correlations are the weakest. Shown in the in- 
set are numerical results for 'pGSE(r) for i?=0.055i?o in 
comparison with the RMT result, Eq. [T51 



V. CONCLUDING REMARKS 

To summarize, we have investigated the statistical 
properties of wave functions in (two-dimensional) chaotic 
nanostructures with spin-orbit interactions, focussing 
particularly on spatial correlations of eigenfunctions. Nu- 
merical results obtained for a chaotic stadium billiard 
were compared with (analytic) results from RMT. It was 
found that excellent agreement between RMT and micro- 
scopic simulations are obtained in a "mean-field" descrip- 
tion of the GOE-GSE crossover. A key observation from 
our results is that correlations of wave function ampli- 
tudes are suppressed with SOC. Interestingly, however, 
these correlations with SOC play a more significant role 
in the two-point distribution function(s). 

Our results have implications for a number of systems 
of current interest. Indeed, the effects of SOC have been 
observed in transport through quantum dots.— These ef- 
fects could also be observed in quantum corrals defined 
on Au and Ag (111) surfaces j2I where large SOC has been 
observed recently^ 2 - especially as scanning tunneling mi- 
croscopy techniques have exquisite control of positioning 
and correlation measurements. 
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